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Abstract. The time it takes a quantum system to complete a tunneling event (which in the 
case of cross-barrier tunneling can be viewed as the time spent in a classically forbidden area) 
is related to the time required for a state to evolve to an orthogonal state, and an observation, 
i.e., a quantum mechanical projection on a particular basis, is required to distinguish one state 
from another. We have performed time-resolved measurements of Landau-Zener tunneling of 
Bose-Einstein condensates in accelerated optical lattices, clearly resolving the steplike time 
dependence of the band populations. The use of different protocols enabled us to access the 
tunneling probability, in two different bases, namely, the adiabatic basis and the diabatic basis. 
The adiabatic basis corresponds to the eigenstates of the lattice, and the diabatic one to the free- 
particle momentum eigenstates. Our findings pave the way towards more quantitative studies 
of the tunneling time for LZ transitions, which are of current interest in the context of optimal 
quantum control and the quantum speed limit. 



1. Introduction 

Tunneling is one of the hallmarks of quantum systems, and physical effects associated with 
quantum tunneling are important in many branches of science [1]. While the probability for 
quantum tunneling can be readily calculated for a variety of systems and has been measured 
experimentally to great accuracy, the time it takes a quantum system to complete a tunneling 
event is a much less well-defined notion. It is often hard to measure in experiment and in many 
cases still the subject of intense debate [2]. In this paper we address the problem of measuring the 
timescale associated with one of the conceptually simplest tunneling phenomena, Landau-Zener 
(LZ) tunneling. 

Landau-Zener tunneling arises when two energy levels of a quantum system cross as a function 
of some parameter that varies in time. There is a possibility of a transition if the degeneracy 
at the level crossing is lifted by a coupling and the system is forced across the resulting avoided 
crossing by varying the parameter that determines the level separation. LZ tunneling was first 
studied theoretically in the early 1930's in the context of atomic scattering processes and spin 
dynamics in time-dependent fields O HI El E] 

In its basic form the LZ problem can be described by a simple two-state model with a 
Hamiltonian given by 
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Figure 1. Energy levels as a function of time. The dashed lines show the so-called diabatic levels i.e., 
the energy position of states in the absence of the interaction. The solid lines demonstrate the so-called 
adiabatic levels, i.e., the eigenstates of the system corresponding to the instantaneous Hamiltonian. 



where the off-diagonal term, A£'/2, is the coupling between the two states and a is the rate of 
change of the energy levels in time. The dynamics of the system can be described either in the 
diabatic or in the adiabatic basis. The diabatic basis is the basis of the bare states of Eq. ([1]) 
when there is coupling (i.e., no off-diagonal entries in the matrix). The adiabatic basis, on the 
other hand, is the basis of a system with a finite coupling AE/2 between the two states. The 

Hamiltonian has two adiabatic energy levels E± = {2atf + /^E"^. 

Assuming that the system is initially, at t — oo, in the ground energy level E- and if 
the sweeping rate is small enough, it will be exponentially likely that the system remains in 
its adiabatic ground state E- at t — ?> -|-oo. The limiting value of the adiabatic LZ survival 
probabilities (for t going from — oo to — oo) is [7J, 



where we have introduced the dimensionless adiabaticity parameter 7 = Aha / /S.E'^ . 

While the above analysis can predict the probability for LZ tunneling very accurately, in 
contains no reference to the dynamics around the avoided crossing and hence to the time it 
takes the system to complete the tunneling event which eventually results in the populations 
measured in the ground and excited bands (or in the corresponding diabatic states). Also, in 
the case of LZ tunneling, which occurs in an abstract space spanned by the energy levels of the 
system as a function of a parameter, the concept of tunneling time is less intuitive than in the case 
of cross-barrier tunneling in real space, to name just one example. Nevertheless, the tunneling 
time (or transition time or jump time, as it is called in certain contexts) associated with LZ is a 
meaningful concept referring to the timescale on which the system evolves around the avoided 
crossing. Analytical estimates for the LZ transition times have been derived in [HI |9] using the 
two-state model of Eq. [TJ In a given basis, e.g., adiabatic or diabatic, different transition times 
are obtained. Vitanov |;9] calculated the time-dependent diabatic/adiabatic survival probability 
at finite times. Analytical estimates for the LZ transition times were derived in [8j using some 
exact and approximate results for the transition probability. 

In general, the LZ jump time in a given basis can be defined as the time after which the 
transition probability reaches its asymptotic value. From this definition one can expect to 
observe a step-like structure, with a finite width, in the time-resolved tunneling probability. 
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Because the step is not very sharp, it is not straightforward to define the initial and final times 
for the transition. It is even less obvious how to define the jump time for both small and for large 
coupling. Some possible choices have been used by Lim and Berry [lOj and Vitanov [8*, '9]. The 
problem is even more complicated when the survival probability shows an oscillatory behavior 
on top of the step structure. The oscillations give rise to an additional time scale in the system, 
namely an oscillation time and a damping time of the oscillations appearing in the transition 
probability after the crossing. Therefore, a measurement of the tunneling time depends very 
much on how these times are defined and also which basis is considered. 
In [S] the jump time in the diabatic/adiabatic bases is defined as 

jump _ -Pd/a(Qo) . . 

where Pd/a. is the transition probability between the two diabatic/adiabatic states, respectively. 
-^d/a(^) denotes the time derivative of the tunneling probability evaluated at the crossing point. 

From this definition, the diabatic iump time is calculated as r^^™"P ~ y^2TTh/a is almost constant 
for large values of the adiabaticity parameter 7. For 7 <C 1, on the other hand, it decreases 
with 7 , rj"™^ ~ 2^/h{^ya)~^. In the adiabatic basis, when 7 is large, the transition probability 
is similar to that in the diabatic basis. For a small adiabaticity parameter, because of the 
oscillations appearing on top of the transition probability step structure, it is not straightforward 
to define the initial and the final time for the transition. Vitanov defines the initial jump time 
as the time t < at which the transition probability is very small (i.e., -Pa(i") = e-Pa(oo), where 
e is a proper small number). The final time of the transition i > is defined as the time at 
which the non-oscillatory part of -Pa('7") is equal to (l + e)Pa(oo). Using these definitions, Vitanov 
derived that the transition time in the adiabatic basis depends exponentially on the adiabaticity 
parameter, rl""^? ^ (i/e)^^^ j-'^/^exp {n/iG-f)) y/W/^. 



2. Experimental results 

The Landau-Zener model is realized in our experiments using Bose-condensed rubidium atoms 
inside an optical lattice I12j . Initially, we created Bose-Einstein condensates of 5 x 10^ 
rubidium-87 atoms inside an optical dipole trap (mean trap frequency around 80 Hz). A 
one-dimensional optical lattice created by two counter-propagating, linearly polarized gaussian 
beams was then superposed on the BEC by ramping up the power in the lattice beams in 
100 ms. The wavelength of the lattice beams was A = 842 nm, leading to a sinusoidal potential 
with lattice constant c^l = A/2 = 421 nm. A small frequency offset Az/(t) between the two 
beams could be introduced through the acousto-optic modulators in the setup, which allowed 
us to accelerate the lattice in a controlled fashion and hence, in the rest-frame of the lattice, to 
subject the atoms to a force Flz = Mclz with olz = ^L^^^^- 

The energy level structure of Bose condensates in optical lattices can be represented by 
energy bands in the Brillouin zone picture. At the edge of the Brillouin zone successive bands 
are separated by gaps, and in the vicinity of the zone edge our system approximates the LZ 
model very well. We can make time-resolved measurements of a single tunneling event in the 
following way: First, the Bose condensate is loaded adiabaticallly into a lattice, after which the 
lattice is accelerated to some finite quasimomentum. Thereafter, the instantaneous populations 
of the eigenstates of the system are measured, the exact protocol depending on the basis chosen. 
In detail, the protocols are as follows: 

• For measurements in the adiabatic basis, after loading the BEC into the optical lattice the 
lattice was accelerated with acceleration clz for a time tLZ- The lattice thus acquired a 
final velocity v = olz^lz- At time t = ti^z the acceleration was abruptly reduced to a 



smaller value Ogep and the lattice depth was increased to V^ep in a time tramp ^ ^B- These 
values were chosen in such a way that at time t = tLZ the probability for LZ tunneling from 
the lowest to the first excited energy band dropped from between ~ 0.1 — 0.9 (depending 
on the initial parameters chosen) to less than ~ 0.01, while the tunneling probability from 
the first excited to the second excited band remained high at about 0.95. This meant that 
at t = ti^z the tunneling process was effectively interrupted and for t > tLZ the measured 
survival probability P{t) = No/Ntot (calculated from the number of atoms A^o in the lowest 
band and the total number of atoms in the condensate A'tot) reflected the instantaneous 
value P{t = ^Lz)- 

The lattice was then further accelerated for a time tsep such that Oseptgep ~ 2nprec/^ (where 
typically n = 2 or 3). In this way, atoms in the lowest band were accelerated to a final 
velocity v ~ 2nprec/Af ; while atoms that had undergone tunneling to the first excited band 
before t = ^lz underwent further tunneling to higher bands with a probability > 0.95 and 
were, therefore, no longer accelerated. At time tgep the lattice and dipole trap beams were 
suddenly switched off and the expanded atomic cloud was imaged after 23 ms. In these time- 
of-flight images the two velocity classes and 2npj-cc/M were well separated and the atom 
numbers A'^o and A^tot could be measured directly. Since the populations were effectively 
"frozen" inside the energy bands of the lattice, which represent the adiabatic eigenstates of 
the total Hamiltonian of the system, this experiment measured the time dependence of the 
LZ survival probability in the adiabatic basis. 

• For measurements in the diabatic basis, after the initial loading phase the lattice was 
accelerated with acceleration clz for a time tLz as in the adiabatic case. At that 
point, however, the atomic sample was projected onto the free-particle diabatic basis by 
instantaneously (within less than 1 /xs) switching off the optical lattice. After a time-of- 
flight the number of atoms in the v = and v = 2pj-cc/M momentum classes are measured 
and from these the survival probability (corresponding to the atoms remaining in the v = 
velocity class relative to the total atom number) is calculated. 

The results of typical measurements in the adiabatic and diabatic bases are shown in Fig. 
[21 The step-like behaviour of the survival probability around t = 0.5Tb is clearly visible, which 
demonstrates that our experimental protocol does, indeed, allow us to access the timescale of the 
LZ transition. It is also obvious from the figure that while in the adiabatic basis the transition 
is smooth and can be well fitted with a sigmoid function, in the diabatic basis there are strong 
oscillations for times t > 0.5 Tg. As a consequence, the tunneling time in the adiabatic basis 
can be easily identified as the width of the transition curve (indicated in the figure), while in 
the diabatic basis it is less obvious when the tunneling event is completed. 

3. Conclusions 

We have demonstrated that experiments with Bose-condensates in accelerated optical lattices 
allow access to the full dynamics of LZ tunneling and hence to the timescales involved in the 
tunneling process, both in the adiabatic and diabatic bases of the problem. Our experiments 
pave the way towards a thorough investigation of tunneling times and the quantum speed limit. 
The latter has recently been discussed theoretically in the context of optimal quantum control 
|13)QSL and more generalized Landau-Zener protocols involving non- linear sweep functions that 
are predicted to lead to shorter minimum times for completing a tunneling event. As our 
experimental setup allows us to realize arbitrary protocols for the lattice acceleration, such 
experiments should be relatively straightforward to realize. 




Figure 2. Timc-rcsolved measurements of LZ tunneling in the adiabatic (filled circles) and diabatic 
bases (open circles)for Fq = 1.197 and y/-Erec = 1-8- The solid line is a sigmoid fit to the adiabatic 
survival probability, and the vertical dashed lines indicate the position of the zone edge at t = 0.5 Tb and 
the tunneling time Ta- 
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